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In a recent post on arXiv [T] Liu & Gunton comment on the role of a bulk chemical potential e in 
determining the radius i? of a self-assembled bundle of chiral protohlaments of indeterminate length 
L ^ R. In essence their concern was that our work ^ erroneously omitted the effect on 72 of a bulk 
chemical potential representing exchange between the solvated pool of monomers and the polymerised 
bundle (polymer). 

We believe that Liu &: Gunton's criticism is fundamentally flawed because they work in the en- 
semble of fixed e and, implicitly, fixed polymer length. In this ensemble e is indeed conjugate to 
R'^ and so can control the polymer radius. However, the physically relevant ensemble is one in which 
the polymer length is free to adjust. Here e is conjugate to L and the aggregate can grow by 
elongating at a fixed equilibrium radius R, corresponding to the minimum polymer free energy density, 
here denoted epoi. When polymers form the monomeric concentration falls and the chemical potential 
adjusts until, finally, e = epoi and only polymers with this equilibrium radius are (marginally) stable. 
Thus, whenever polymers form, the final value of e is not a free parameter. There is a simple 
physical argument for why the radius of the aggregate must be insensitive to its length, and hence to 
the total mass of aggregated polymer: In the absence of long ranged interactions there is no physical 
coupling between R and L. 

Since our original letter other authors have undertaken similar studies in which a bulk chemical 
potential was also included [3]. These authors take the limit in which the polymer length L —?■ oo, 
eliminating what is otherwise an apparent dependence of R on e, e.g. in Eq 7 in 

In the interests of completeness we would add that the chemical potential can affect the aggregate 
radius but only if the concentration is either near the CMC, where the aggregates are so small as to 
no longer be essentially linear, or so high that strong polymer-polymer interactions are important. 
Neither regime is the focus of [1], [2] or [3]. 

For those with an interest in the total polymer length, equivalent to the polymer concentration, we 
add that this could be calculated from epoi given in [2] (or [3]), combined with an appropriate relation- 
ship between e and the concentration of monomers in solution (j). The final monomer concentration is 
given by e((/>) = Cpoi and the corresponding increase in the amount of polymer follows by conservation 
of material. 
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